DETERMINISTIC RATCHETS: ROUTE TO DIFFUSIVE TRANSPORT 
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The rectification efficiency of an underdamped ratchet operated in the adiabatic regime increases 
according to a scaling current-amplitude curve as the damping constant approaches a critical 
threshold; below threshold the rectified signal becomes extremely irregular and eventually its time 
average drops to zero. Periodic (locked) and diffusive (fully chaotic) trajectories coexist on fine 
tuning the amplitude of the input signal. The transition from regular to chaotic transport in 
noiseless ratchets is studied numerically. 
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I. INTRODUCTION 

Particles in an asymmetric potential can drift on aver- 
age in one direction even when the time and space aver- 
age of all applied macroscopic forces or gradients is zero. 
To achieve directed transport in such a device, called a 
ratchet, a time-correlated source of energy is required, 
for instance non-markovian external fluctuations (ther- 
mal ratchets 0) or a time periodic drive (rocked ratch- 
ets ^]). In signal analysis notation, a ratchet can be re- 
garded as a (random or periodic) signal rectifier. Rocked 
ratchets with massive particles exhibit strong inertial ef- 
fects capable of reversing their current ||. Moreover, in- 
ertial rocked ratchets are naturally subject to developing 
chaotic dynamics, the onset of which becomes detectable 
as thermal fluctuations are switched off [0, || . 

This class of devices, operated at zero noise level, may 
be assimilated to an asymmetric version of the damped- 
driven pendulum, a chaotic system investigated at depth 
in the Eighties M-S- The dynamics of a massive par- 
ticle in a cosine potential was reproduced in terms of a 
"climbing-sine map" 0: Running orbits, relevant to the 
rectification mechanism in a ratchet, can be either pe- 
riodic or diffusive, depending on the value of the map 
control parameter (viz. the amplitude of the sine term). 
The phase-space portrait of the actual damped-driven 
pendulum was computed by Huberman et al. ||, who 
revealed the existence of delocalized strange attractors 
with an intricate structure on all scale, later recognized 
to be fractal objects ||. 

In this article we investigate the transport of unit mass 
particles in a ratchet potential subjected to a sinusoidal 
driving force A(t) with angular frequency much smaller 
than the damping constant 7 and the librational fre- 
quency luq of a particle trapped at the bottom of a po- 
tential well. The damping constant is lowered from ex- 
ceedingly large values 7 > IOcjo, where the Smoluchowski 
approximation of Rcf. |J holds good, deep into the un- 
derdamped regime 7 ~ O.Olwo. The adiabatic condition 
£1 < c<Jo,7 imposed here, has not been addressed in the 



damped-driven pendulum literature, but is believed to 
maximize transport in biological systems jj] and Joseph- 
son junction arrays |i§ ||, [n| , among others. On decreas- 
ing 7 we noticed that (a) The rectification efficiency of 
an inertial ratchet increases inversely proportional to 7 
according to a scaling current-amplitude (v)-Aq curve; 
(b) Below a critical 7 value, more (v)-Aq curves coexist, 
thus signaling the appearance of regular trajectories with 
higher periodicity; (c) On further decreasing 7, all trajec- 
tories develop a finite diffusion coefficient and the ratchet 
current drops monotonically toward zero. These results 
are of immediate application to the design and operation 
of nano-array devices, both in solid state- jjj, |l0| and in 
bio-engineering JjJ. 



II. THE DOUBLE-SINE MODEL 



mo 



We have integrated numerically the archetypal ratchet 
del § 



x = -jx - V'(x) + A sin(nt), 
with asymmetric double-sine potential (Fig. 



V(x) 



[sin(fcx) H — sin(2fcx)] 



(1) 



(2) 



and arbitrary time origin of the drive A(t) — Aq sin(fit). 
The constant k was set equal to 7r/2, so that the unit cell 
of (0) is I — 4 and the libration frequency at the bottom 



of its wells is w = \J fc(3' V / 3/2) 1 /2 ~ 1.591. The recti- 
fication window (Ai,A2) of the ratchet (Q)-(|2j) is deter- 
mined by the depinning thresholds of the tilted potentials 
V(x)tx\A \, i.e. Ax = 3/4 to the right (-) and A 2 = 3/2 
to the left (+). As A\ < A2, the ratchet current is nat- 
urally directed to the right Q or, equivalently, the time 
average (v) of the particle velocity is positive definite for 
Aq > A\, The forcing frequency f2 is taken so small that 
the the fi-dependence of A±, A2 is negligible H. 
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In the adiabatic limit, £1 <§; a»o,7, the values of (v) 
are quantized with step Av — I/Tq, corresponding to 
a net particle shift of one unit cell I per forcing period 
Tq = 2it/il. Of course, quantization of the (v)-Aq char- 
acteristics becomes less and less apparent as the ratchet 
operating conditions most relevant to the present inves- 
tigation, i.e. fi <C 7 <C u)q, are approached (see Fig. ||). 
The integration of Eq. (|l]) was performed by means of a 
4-th order Runge-Kutta algorithm with integration step 
At = 2- 10- 6 T n . 

In our numerical work we set Q = 0.01 and explored 
the ratchet response on varying the damping constant 
in the interval [10 _2 ,10]. For large 7 values (7 = 10 
in Fig. ||) we reproduced the results of Ref. 0. The 
current-amplitude curve grows linearly with A from zero 
up to a maximum (A 2 — Ai) 2 /2, between A\ and A2, and 
then decays with comb-like structure reminiscent of the 
Vernier effect. [The amplitude of the particle excursions 
both to the left and to the right increases fast with Aq, 
but slightly faster to the left than to the right.] Note 
that the rectification effect persists for values of Aq much 
larger than A 2 . 



A. The Mobility curves 

For 7 values in the range [0.07, 0.5], the (v)-A charac- 
teristics exhibits the universal behavior displayed in Fig. 
||. The data points of 7(1;) versus Ao rest on a unique 
limiting curve (7 = 0.1 in Fig. |^) with deviations of a few 
percent (and diminishing with lowering 7). The scaling 
law (v) oc I/7 for an adiabatically driven underdamped 
particle is to be expected Jll], [l^] . More remarkably, we 
observed that: 

(i) The ratchet current gets more and more narrowly 
restricted to the rectification window (Ai,A 2 ) as the 
damping constant grows smaller than luq. This effect is 
even more apparent in Fig. || for 7 = 0.05; 

(ii) The onset of the ratchet current at the threshold 
A\ is abrupt; 7(1*) jumps from zero up to 0.086 within one 
incremental step AAq = 2 ■ 10~ 6 . On the contrary, the 
sharply decaying branch of the curve "/(v) for Aq > A 2 is 
continuous, its persistent irregularity being a signature 
of the underlying chaotic dynamics of the system (see 
caption of Fig. ||); peaks and dips have been resolved on 
the scale AAq = 10~ 5 , i.e. within the thickness of the 
broken line connecting our data points. 

The universal j(v)-Aq curve can be determined as fol- 
lows. For A\ < Aq < A 2 the particle runs to the right 
with velocity v(t) = Ao sin(f2< +(/))/ where the adiabatic 
phase-lag is <j) — — arctan(f2/7). As discussed in Ref. 
JTi"[ , a trapped particle gets depinned to the right when 
A(t) > A\ and, then, repinned as soon as A(t) < A r , with 
A r oc Ai(7/wo) — 0. An underdamped particle will be 
then traveling in the positive direction during the time in- 
terval [fi -1 arcsin(Ai/Ao), 7r/f2] mod(27r/J7); hence, the 



time average of v(t) over one forcing period yields 



ee »(A , Ax) = —[y/i-^/AoF + l]. (3) 

Aq ZTT 

Ther.h.s. ofEq. (||) is 7 independent and accurate in the 
limit 7/wo - * 0. For Aq > A 2 a similar argument leads 
to j{v)/A = fi{A ,Ai) - fi{A ,A 2 ) (see Fig. §, inset). 
The analogy with the experimental I — V characteristics 
of the Josephson ratchet in Ref. (HJ is apparent. 

A totally different scenario emerges on decreasing 7 
below 0.07. The current for A > A 2 drops close to zero, 
while in the rectification window the data points taken at 
steps of AAq = 10~ 3 arrange themselves along different 
current-amplitude characteristics. In particular: 

(iii) In Fig. |^, besides the universal curve (||) denoted 
by 1 : 1, at least three more characteristics, 1:2 — 1:4, 
are resolved. In our notation curve 1 : m corresponds to 
curve (||) divided by the integer m. Note that, contrary 
to the case of the damped-driven pendulum j^, ^ , in the 
asymmetric problem (|l|)-(^) odd m values are not ruled 
out; 

(iv) A j(v)-A background is enclosed by curve 1 : 4. 
(or, possibly, by 1 : 5, not shown) The data point distri- 
butions on the characteristics 1:1 — 1:4 are denser in 
the vicinity of A\ , while the background gets thicker close 
to A 2 . [Sparse negative mobility background values have 
been recorded immediately above the threshold A\\ In 
the window (j4i, A2) compact Aq intervals corresponding 
to the diverse characteristics (or to the background) can 
be easily indentified on the scale A^4o = 10~ 5 (see Fig. 
||, inset). On further decreasing 7 below 0.04, the curves 
1:1 — 1:4 get depleted and only an extremely irregular, 
low background survives. 

In order to explain the coexistence of many (v)-Aq 
characteristics, we sampled one trajectory from each of 
the uppermost three 1 : m curves (Fig. top panel, 
upper). These trajectories are seemingly periodic and 
independent of the initial conditions. The trajectory cor- 
responding to curve 1 : 1 drifts with positive speed step- 
wise, by alternating one plateau (trapped particle) with 
one jump of constant length Ax p (running particle). The 
trajectory representative of curve 1 : 2 consists of a reg- 
ular sequence where each plateau is followed by a full 
oscillation (here, with measured peak-to-peak amplitude 
L ~ 3800). Note that the separation Ax p between two 
adjacent plateaus is almost the same as in the first case, 
being the Aq difference as small as 10~ 4 . Consequently, 
the average velocity of a type 1 : 2 trajectory is half the 
reference value (||) on curve 1:1. Analogously, a type 
1 : 3 trajectory alternates periodically one plateau each 
two oscillations, thus drifting with one third of the veloc- 

ity (§■ 

B. The Diffusion Coefficient 

The behavior of the background trajectories is even 
more interesting: 
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(v) Plateaus are separated by random sequences of full 
oscillations, extremely sensitive to the initial conditions 
(Fig. f|, top panel, lower); the lengths r of such sequences 
are distributed according to a simple exponentially decay- 
ing law. Occasionally, the oscillation sequences slant to 
the left, i.e. in the negative direction, with glitches of dif- 
ferent length. As suggested in Re fs. |(| , we determined 
the second cumulant of the random trajectories plotted 
in Fig. |, 

C(t) = x^) -W) 2 , (4) 



the averages (...) being taken over the particle initial 
conditions x(0),x(0) 0. From the fitting law 

lim C(t) = 2Dt, (5) 

t — >oc 

one eventually extracts the diffusion coefficient D of the 
noiseless process x{t) (Fig. EL bottom panel); 

(vi) Trajectories change from regular to chaotic on tun- 
ing 7 at constant A®. Typically, for values of Aq corre- 
sponding to points sitting on curves 1 : 2 — 1 : 3 in Fig. 
|3|, the periodicity of the x(t) trajectories varies as 7 is 
ramped from 0.07 down to 0.04; trajectories with dif- 
ferent periodicity show up repeatedly, interspersed with 
diffusive trajectories; eventually, the diffusive nature of 
the underdamped processes prevails. Reentrant type 1 : 1 
trajectories do not imply the current reversals observed 
in Ref. g. For £1 = 0.01, type 1 : to trajectories with 
to > 4 are hardly detectable; moreover, these different 
behaviors turned out to be insensitive to the initial con- 
ditions. 

The onset of diffusive transport is a main feature of 
inertial ratchets Although an analytical estimate of 
the diffusion coefficient D lies beyond the grasp of our 
theoretical tools a simple relationship between drift and 
diffusion of a ratchet (|l|)-(§) in the fully chaotic regime 
can be easily established. Consider the case A) = 0.9305 
represented in Figs. || and [J. The average value of j(v) 
is 0.045, namely one fifth of the corresponding value on 
curve 1 : 1; thus, the average duration of the oscillation 
sequences along the relevant diffusive trajectories is of 
the order f = uTq with n — 4. The separation Ax p 
between two subsequent plateaus [no matter what n, as 
long as one disregards the negative slant of the oscilla- 
tion sequences] can be evaluated by integrating v{t) — 
A sin(rtt + along the running solution branch, i.e. 
over the time interval [f2 _1 arcsin(A 1 /^4 ), for the 

parameter values of Fig. ||, Axj, ~ (3/4)L. As a conse- 
quence, the standard formula [ |l3f 

D = (Ax p ) 2 ^-=f^ (6) 

yields the estimate D = (3L/4) 2 /12r ~ 270 that com- 
pares fairly closely with the fitted value D = 380 ± 50 (in 
the units of Fig. [|) . 



III. PHASE-SPACE PORTRAITS 

The chaotic nature of the trajectories in Fig. [| is il- 
lustrated by the reduced phase-space portrait of Fig. ||; 
here, 7W and x = xmod(Z) — (1/2) have been recorded 
once every forcing cycle at t n = (n + l/4)Th (strobo- 
scope snapshots). The periodicity of type 1:1 — 1:3 
trajectories is apparent: In all three cases there exists 
one single point at a lowest velocity vq , associated to the 
depinning of the ratcheted particle out of the potential 
trap. Type 1 : 2 and 1 : 3 trajectories exhibit one or two 
full oscillations, respectively, prior to repinning (Fig. [J); 
therefore, the snapshot velocities corresponding to such 
oscillations are larger than the relevant vq. In view of a 
straightforward energetic argument, one concludes that 
the portrait points representative of regular trajectory 
oscillations must be located along the upper asymmetric 
band of Fig. ||, v\ (x) , defined by 

Vi(x) - vi= -V(x-xfi)/vj, (7) 

where i = 1, 2 is a band index, u[ its mean value, V(x) the 
ratchet potential (|2|) and xq a spatial shift that depends 
on the definition of x and on the velocity phase-lag 4>. 

Most remarkably, the phase-space portrait represen- 
tative of regular trajectories consists of singular points 
(more likely, confined segments), denoted by 1-3 in Fig. 
H this implies that such trajectories are truly periodic 
and that, equivalently, for appropriate parameters values 
inertial depinning amounts to a renewal process JD| of 
variable locking time. The velocity gap v\ — vq in our stro- 
boscope rappresentation can be easily estimated. Let us 
assume that, on neglecting a small phase-lag 0, the par- 
ticle velocity is zero at t = (no tilt) and relaxes towards 
its asymptotic value A(ti)/"/ — Aq/j according to an ex- 
ponential law, namely, vq = (A /^)[l — e -7 ^' 1- * ^], if the 
particle exits a potential well assisted by the tilt A(t) 
at t = fl- 1 a,TCsm(Ai/A ), or W[ = (A /7)[l - e~ 7tl ], 
if, thanks to inertia, it executes a full oscillation with- 
out getting trapped. The ensuing estimate 7(1^"— ^o) — 
3.8 • 10~ 2 agrees fairly well with our data in Fig. ||. 

The reduced phase-space portrait of a diffusive trajec- 
tory is characterized by a confined segment ("a" in Fig. 
||), associated with the depinning events, and a compli- 
cated gapped band structure, accounting for the negative 
trajectory slant. The splitting of the upper edge into two 
close bands Vi (x) with i = 1 , 2 is related to the existence 
of two flexural points on the shallow, rising branches of 
V(x). The dispersion of the representative points along 
the bands Vi(x) and the bottom ribbons (same shape as 
Vi(x), but apparently quantized along the 7W axis) is a 
clear-cut signature of a fully chaotic behavior. 



IV. CONCLUSIONS 

In this report we focused on certain properties of adi- 
abatic inertial ratchets that are of direct application to 
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bio-engineered and nanoelectronic devices. Related ques- 
tions that might deserve closer attention by theorists and 
experimenters, alike, are: (a) The robustness (or selec- 
tivity) of curves 1 : m in the presence of noise. In par- 
ticular, one wonders if all characteristics 1 : m may sur- 
vive in the presence of thermal fluctuations. Most likely, 
with increasing the temperature, the diffusive nature of 
the ratchet trajectories will prevail in the vicinity on A2, 
first, and eventually over the entire rectification window, 
(b) The dependence of the diffusion coefficient on the 
ratchet parameters. In the absence of noise, the diffusion 
coefficient depends crucially on the nature of the rele- 



vant trajectories. Trajectories may be periodic (D = 0), 
fully chaotic (when the argument of Eq. @ applies) or 
exhibit an intermediate behavior, termed confined chaos. 
For trajectories in this category all singular points in the 
phase-space portraits of Fig. H are replaced by curves of 
different length and the diffusion coefficient acquires a fi- 
nite value, (c) The generalization of the present analysis 
to wider classes of asymmetric potentials [|l2[ |f4) . While 
the mobility curves may depend on the choice of V(x), 
the dynamical regime Q <C 7 <C 1^0 is expected to show 
important similarities, no matter what the details of the 
substrate model. 
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FIG. 1: The double-sine potential (Q) with k = ty/2. As the l.h.s. of each potential well is steeper than the r.h.s., the natural 
ratchet direction is positive. 
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FIG. 2: j(v)-Ao characteristics of ratchet (|l|)-(^|) for Q, — 0.01 and different values of 7 (incremental steps AAo = 10 -3 ). 
The dashed curve represents the adiabatic limit = 10~ 4 and 7 = 10. Inset: the numerical data for 7 = 0.1 (solid curve) are 
compared with fi(Ao,A\) for Ai < Ao < A2, Eq. (^|), and fJ,(Ao, Ai) — ^{Aq, A2) for Ao > A2 (dashed curve). The velocity 
averages have been taken over 300 forcing cycles, after discarding the first 10 3 cycles to get rid of transient effects. 




FIG. 3: "/{v)-Ao characteristics of ratchet (M-fg) for Q = 0.01, 7 = 0.05 and AA = 10" 3 (circles). Inset: details with 
AA = 10~ 4 . The curve 7 = 0.1 from Fig. 13 is plotted for reader's convenience (solid curve). Averages have been taken 
as in Fig. [2], but for the squares on curve 7 = 0.1, where the preparation time was doubled (2 • 10 3 cycles). Dashed curves: 
characteristics 1:1-1: 4. 




1 2 3 4 5 

FIG. 4: Top panels: samples of regular and diffusive trajectories (the latter ones for different initial conditions, see text) for 

Q. = 0.01 and 7 = 0.05; bottom panel: C(t) — x 2 (t)~x(t) for Aq — 0.9305; here, (. . .) are taken over 10 2 different preparations. 
Dashed line: fitting law (^J) with D = 380. 
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FIG. 5: Reduced phase-space portrait for A = 0.9309 (1), 0.9308 (2), 0.9307 (3), and 0.9305 (dots). Dashed straight lines 
connect the portrait point sets of (2) and (3). Open circles: vi(x) (upper) and V2(x) (lower), Eq. (Q), with xo = —0.11. Note 
that the dots in the diffusive regime are concentrated at the confined segment "a" and on the two upper bands. The remaining 
parameters are as in Fig. ^[ 



